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^O \ Abstract. Planetary nebulae (PNe) are very important kinematical 

tracers of the outer regions of early- type galaxies, where the integrated 
light techniques fail. Under ad hoc assumptions, they allow measurements 

QQ ' of rotation velocity and velocity dispersion profile from descrete radial ve- 

f^ , locity fields. We present the results on the precision allowed by different 

set of radial velocity samples, discuss the hypotheses in the analysis of 
descrete velocity fields and their impact on the inferred kinematics of the 

Q^ ' stellar population. 

' ' 1. Introduction 

"^ I With 4m telescopes and multi-object and/or multifiber spectrographs it has been 

Ci • possible to obtain information on the kinematics of the outer regions of ellipticals 

by measuring the radial velocities of individual stars during their phase of plane- 

- tary nebulae. This is crucial because PNe are found at very large distances from 

K> \ galaxy centers where kinematical measurements based on standard integrated 

j^ ■ stellar light techniques are no longer possible. These radial velocity samples 

obtained from the PNe 5007 A[OIII] emission in giant Es and SOs (D=15-17 

Mpc) contain up to 50 PNe radial velocities (Arnaboldi et al. 1998), and larger 

samples are available only for nearer objects (i.e. NGC 5128, Sombrero, NGC 

3115). Therefore there might be biases introduced by small number statistics 

that need to be investigated and understood. 

Based on a statistical approach, we investigate the possibility of re-building 
the actual kinematics of spherical systems starting from descrete radial velocity 
fields. We build equilibrium systems for which dynamical and kinematical pa- 
rameters are known (spherical model -|- known velocity field). By Montecarlo 
simulations we produce 100 "observational sets", each for a given sample size, 
i.e. 50, 150 or 500 randomly chosen stars, and then analyse each sample to de- 
termine the rotation curves and velocity dispersion profiles. We account also for 
the measuring errors (~30% of the maximum rotational velocity). The aim of 
this work is to address the following questions: what is it the precision allowed 
by the different statistical samples in determining the kinematical quantities? 
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Furthermore, may the hypotheses on the internal rotational structure introduce 
any biases? 



2. Model procedure 

We assume a constant M/L ratio and the Hernquist model (Hernquist 1990) for 
the dimensionless mass density distribution 

where r is in units of a core radius a and A^ is a normalisation constant chosen 
as to have a total mass M(=l within a distance of 18a from the center (~ 10i?e)- 
The cumulative mass distribution is 

Mi{x) = An r x'^p{x')dx' (2) 

Jo 

In spherical models with dark matter, we consider an additional mass contribu- 
tion coming from an isothermal halo, and can write the dark matter density as 
(Dubinsky & Carlberg 1991): 

Mgd 1 
Pd{r) = , , ,.3 (3 

zvr r (r + dy 

where M^, d and r have the same meaning as in eq. (|l|). In eq. (^), we take 
d = 10a and AI^ = 1.1. Mi. The cumulative mass distribution is defined as in 
eq. (^) with total mass: 

Mt{r)=Mi{r)+Md{r). (4) 

The model for the rotational velocity is 

Vrotix) = 



r'n + x^ 





where Vb is the maximum rotational velocity, ro is a scale factor and x is the 
distance from the rotational axis (cylindrical rotation). 

The PNe distribution is computed taking into account the selection effects of the 
bright continuum background from the stars in the central region of the galaxy. 
This selects a projected radius Kum onto the sky plane out of which the PNe 
sample is complete. 

3. The descrete radial velocity field 

From the mass density in space, p(r), we extract via Montecarlo a star at posi- 
tion {xp,yp,Zp) = {rp,9p,(pp). We check if Rp = Xp + Yp, where {Xp,Yp) are 
the star projected coordinate onto the Sky Plane, is greater that the Rum, the 
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Figure 1. A simulated radial velocity field: the dotted line indicates 
the direction of niaxiniuni gradient. 



completeness limit radius. In the 3D space, we assign to the selected star its ve- 
locity vector under the hypothesis of ^Hocal isothermal approximation" , i.e. the 
observed distribution of radial velocity profiles in Es are Gaussian with max- 
imum deviation of 10% (Bender et al. 1994). This is consistent with taking 
F(r,v) as product of three Gaussians. The O"r,cre,cr0 components are derived 
from the Jeans equation with the adopted p(r) for the mass distribution, assum- 
ing isotropy for the velocity ellipsoid and the adopted rotation curve. 
Once the 3D velocity vector is computed, the radial velocity is derived by pro- 
jection of the velocity vector along the line-of-sight. This procedure is iterated 
for all the stars in a given sample and the 2-D descrete velocity field VobsiX, Y) 
is derived; a 500 PNe sample is shown in fig.l. 



4. The analysis procedure 

For each given Yobsi^^Y) sample, we perform a simple analysis of the velocity 
field as follows: we fit a rigid rotator (bilinear function) Vrad = « + bX + cY and 
the results of this fit are compared with the results of a flat rotational curve fit 
(flat-fit) of the form Vrad = Vgys + vq cos((^ — c^^o) then we obtain the field of the 
residuals for each of the two interpolated fields Aw = Vobs — Vrad- 
The aim of this procedure is to determine the characteristics of the velocity field: 
systemic velocity, direction of the maximum velocity gradient (Zl), velocity gra- 
dient and maximum rotational velocity. Along Zl, we select a slice of width 
dZl in order to trace the kinematics along this relevant axis; dZl must be small 
enough to trace the kinematics related to Zl and large enough to allow a signif- 
icant statistical sample. To describe the velocity dispersion profile, we analyse 
the residual fields from the two adopted fits, by computing averages along Zl in 
bins, which contains at least 10 points. 

As a check for the presence of biases, along the direction Zl of maximum gra- 
dient, we obtain the rotation curve and the velocity dispersion profile simply 
as the average and radial velocity RMS in each of the given bins used in the 
analysis of the residuals, and we refer to this as the "No-fit" procedure. 
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Figure 2. The grey points are the values of the difference between 
the velocity dispersion determined in the last bin along Zl and the 
expected value (weight average). The continuous line indicates the 
linear regression. The behaviour showed by the residuals indicates 
that higher rotation in the model causes the bilinear fit to introduce a 
larger bias in the velocity dispersion determination. 



5. Results and conclusions 

As a first result we obtained that the direction of maximum gradient is inde- 
pendent of the rotational structure adopted for fit. For simulated samples with 
500 PNe, the comparison between the simulated results with the expected model 
values show that the bilinear fit may introduce an over estimate of the velocity 
dispersion in the outer bins, depending on the intrinsic rotational structure 
of the galaxy. This bias is a function of the v/a ratio as shown in fig. 2. The 
no-fit procedure does not introduce biases in the determination of the velocity 
dispersion profile. For samples with 50 PNe, in the outer bins we over-estimate 
the velocity dispersion by an amount which is dependent from the underlying 
mass distribution. If we take into account the weight effects in the binning, 
using the surface brightness distribution (i.e. surface density distribution via 
M/L ratio), we found a better consistency with the simulated data. The preci- 
sions obtained for the kinematical quantities depend on the selection performed 
along Zl. Our best mean estimate is 15% for 500 PNe and 22% for 50 PNe. 
PNe can be used very efficiently as kinematical probes to trace the dynamics 
of the outer stellar halos in giant early-type galaxies provided that the analysis 
of the descrete radial velocity fields avoids any systematics depending on the 
assumption of the internal angular momentum distribution. 

References 



Arnaboldi, M., Capaccioh, M., Freeman, K.C. et al. 1998, ApJ, 507, 759 
Bender, R. Saglia, R.P., Gerhard, O.K. 1994, MNRAS, 269, 785 
Dubinsky, J., Carlberg, R.G. 1991, ApJ378, 496 
Hernquist, L., 1990, ApJ, 356, 359 



